Single-molecule stochastic resonance 
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Stochastic resonance (SR) is a well known phenomenon in dynamical systems. It consists of 
the amplification and optimization of the response of a system assisted by stochastic noise. Here 
we carry out the first experimental study of SR in single DNA hairpins which exhibit cooperatively 
folding/unfolding transitions under the action of an applied oscillating mechanical force with optical 
tweezers. By varying the frequency of the force oscillation, we investigated the folding/unfolding 
kinetics of DNA hairpins in a periodically driven bistable free-energy potential. We measured several 
SR quantifiers under varied conditions of the experimental setup such as trap stiffness and length 
of the molecular handles used for single-molecule manipulation. We find that the signal-to-noise 
ratio (SNR) of the spectral density of measured fluctuations in molecular extension of the DNA 
hairpins is a good quantifier of the SR. The frequency dependence of the SNR exhibits a peak at 
a frequency value given by the resonance matching condition. Finally, we carried out experiments 
in short hairpins that show how SR might be useful to enhance the detection of conformational 
molecular transitions of low SNR. 

PACS numbers: 



I. INTRODUCTION 



All nonlinear systems that exhibit stochastic noise are 
susceptible to undergo stochastic resonance (SR). When 
SR is triggered, the response of a system to an external 
forcing is amplified. SR has been studied in a large vari- 
ety of systems, including climate dynamics 0,0, colloidal 
particles jj-Hl, biological systems 0-131, and quantum sys- 
tems [1, . With the recent advent of single- molecule 
techniques, it is nowadays possible to measure SR at the 
level of individual molecules. Biomolecules exhibit rough 
and complex free energy landscapes that determine fold- 
ing kinetics and influence the way they fold into their 
native structures. The use of force spectroscopy tech- 
niques has become important practice in studies of molec- 
ular biophysics. By applying a mechanical force at both 
extremities of an individual molecule and by recording 
the time evolution of the molecular extension (the reac- 
tion coordinate in these experiments), information about 
the folding reaction can be obtained. The application 
of forces makes possible to disrupt the weak bonds that 
hold their native structure to reach a stretched unfolded 
conformation. In this way thermodynamics (e.g. the free 
energy of folding) and kinetics (the rates of unfolding and 
folding) can be determined. 

Although most SR studies use temperature as a tun- 
able parameter, this is not the best choice to investigate 
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SR effects at the single-molecule level. Biomolecules have 
a strong sensitivity to temperature variations. Indeed, 
beyond increasing thermally assisted noise, temperature 
also modifies the shape of the molecular free energy land- 
scape. Thus, another tunable parameter such as the os- 
cillation frequency of force might be more appropriate 
to study SR in biomolecules. SR appears as a maxi- 
mum in the response of a biomolecule at a character- 
istic frequency (the resonance frequency). This occurs 
when a characteristic timescale of the signal (e.g. its 
decorrelation or relaxation time) matches half period of 
the oscillation (the so-called matching condition). The 
matching condition must not be taken as a strict equality 
but a qualitative relationship between the two timescalcs 
[ill [T2 |. This means that different SR quantifiers may 
not give coincident resonance frequencies specially for low 
quality resonance peaks. It seems important to investi- 
gate which SR quantifier is best suited to identify SR 
behavior. 

In this work, we use optical tweezers to investigate SR 
in single DNA hairpins driven by oscillatory mechanical 
forces. The high chemical stability of DNA makes DNA 
hairpins excellent models to investigate SR at the single- 
molecule level. When force oscillates around the average 
unfolding force, thermally activated hopping kinetics be- 
tween the folded (F) and unfolded (U) states synchronizes 
with the frequency of the external driving force, leading 
to SR. SR can be measured by recording the oscillations 
produced in the molecular extension, relative to the mag- 
nitude of the noise produced by the thermal forces. Our 
aim in this work is to perform a systematic study of SR 
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FIG. 1: Experimental setup, hopping and SR experiments, (a) Illustration of the experimental system (left), DNA 
sequence of the two-state hairpin HI (upper right, sequence shown in color code) and experimental force-distance curve for HI 
obtained from a pulling experiment (lower right), (b) Typical force and extension traces of the hopping experiments for HI 
obtained in CFM (upper) and PM (lower), (c) Probability distributions of the residence times for HI in the F (red) and U 
(blue) states obtained from the hopping experiments at / = /c ~ 14.5 pN in the CFM. The black curves show the exponential 
fit, (1/a) exp(— r/a), to the data, with a=1.42 s (a=1.34 s) for the F (U) state, (d) Typical force and extension traces (upper 
and middle) obtained by applying a force oscillation protocol with amplitude A = 0.7 pN and frequency Vos = 0.4 Hz around 
the coexistence force /c ~ 14.5 pN. In the lower panel, we show the measured power spectrum, S(i'), calculated by Fast Fourier 
Transform with window size A'^ — 2^^ of X{t) shown in the middle panel. The sampling rate of the instrument is 1 kHz. The 
red area is the output signal (OS, Eq[2]) and the blue vertical bar represents the background noise (BN, EqlSjl. 



in single-molecules exhibiting bistable dynamics, rather 
than using SR as a useful tool to determine the kinetic 
properties of DNA hairpins. In fact, these can be esti- 
mated by using other much less time-consuming methods 
(e.g. by directly analyzing hopping traces). Yet, we also 
carry out SR studies in short hairpins that show how SR 
might prove useful to enhance the detection of conforma- 
tional transitions of low SNR. 

The paper is organized as follows. In Section II, our 
experimental set up is explained. Our main SR results in 
DNA hairpins are presented in Section III, and the influ- 
ence of the experimental conditions (i.e. dsDNA handle 
length and trap stiffness) is investigated in Section IV. 
We compare different SR quantifiers in Section V and in 
section VI we describe the related phenomenon of res- 
onant activation. Finally in Section VII, we purposely 
designed short DNA sequences to increase the noise of 
the signal to test whether SR can still be used to identify 
the hopping rate. In the last section, we summarize our 
conclusions, and discuss situations where SR might be a 



useful technique. 

II. EXPERIMENTAL SETUP AND HOPPING 
EXPERIMENTS 

In Fig. la, we show a schematic illustration of our ex- 
perimental setup (left) and the DNA sequence of hairpin 
HI that we investigated (upper right). The DNA hairpin 
is tethered between two short dsDNA handles (29 bp) 
that are linked to micron-size beads [l^. One bead is 
captured in the optical trap whereas the other is immo- 
bilized at the tip of a glass pipette [l^ . By moving the 
position of the optical trap relative to the pipette, a force 
is exerted at the extremities of the hairpin. 

In a pulling experiment, the optical trap is moved away 
from the pipette and mechanical force is applied to the 
ends of the DNA construct (DNA hairpin plus DNA han- 
dles) until the value of the force at which the hairpin 
unfolds is reached. In the reverse process, the trap ap- 
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preaches the pipette and the force is relaxed until the 
hairpin refolds. In this experiment, the force exerted 
upon the system, /, is recorded as a function of the 
relative trap-pipette distance giving the so-called force- 
distance curve (Fig. la, lower right). Around the co- 
existence force, /c ~ 14.5 pN, the hairpin hops between 
the F and U states for sufficiently low pulling speeds. 



Hopping experiments can be done in two different 
modes: constant force mode (CFM) and passive mode 
(PM) [H [H. In the CFM, the force applied to the 
DNA construct is maintained at a preset value by mov- 
ing the optical trap through force-feedback control (Fig. 
lb, upper). The folding and unfolding transitions of the 
DNA hairpin are followed by recording the trap position, 
X{t). In contrast to the CFM, the PM is operated by 
leaving the position of the optical trap stationary with- 
out any feedback. The bead passively moves in the trap 
in response to changes in the extension of the DNA con- 
struct (Fig. lb, lower). When the hairpin unfolds, the 
trapped bead moves toward the trap center and the force 
decreases; when the hairpin folds, the trapped bead is 
pulled away from the trap center and the force increases. 
The folding and unfolding transitions of the DNA hair- 
pin are followed by recording the force, f{t). In both 
cases (CFM and PM), the kinetic rates of hopping can 
be measured from the residence times of the trace {X{t) 
in the CFM and f{t) in the PM). Fig. lb shows hopping 
traces measured in the CFM and PM at the co-existence 
force, fc — 14.5 pN, where the hairpin hops between the 
F and U states populating them with equal probability 
(i.e. it spends equal time in both states). 



In this work, we focused on the experiments at con- 
trolled force, rather than at fixed trap position. Both 
the hopping and the oscillation experiments (described 
below) were carried out using the force feedback con- 
trol. The reason is that the controlled force experiments 
avoid undesirable drift effects in force that strongly af- 
fect the kinetics of the hairpin (see Methods). Therefore 
we mainly carried out the experiments in the CFM by 
recording the position of the trap, X{t). This signal ex- 
hibits dichotomic motion between the two distinct levels 
of extension (Fig. lb, upper left). The difference between 
the two levels (short extension, folded; long extension, 
unfolded) reflects the release in extension (~ 18 nm) of 
the 44 nucleotides of hairpin HI. From X{t) we can ex- 
tract the residence time distribution at each state that 
shows the exponential form characteristic of first-order 
decay processes (Fig. Ic). The fit of the time distribu- 
tion to an exponential function allow us to get the average 
residence time. The force-dependent kinetic rates (equal 
to the inverse of the mean lifetimes), fcpu a-nd /cuf, were 
measured at the co-existence force, fc — 14.5 ± 0.3 pN, 
giving fcc = fc^u = ^UF 
SI). 



III. SR EXPERIMENTS 

To induce the SR phenomenon, we applied an oscil- 
lating force, f{t), to the DNA hairpin using the force 
feedback protocol, where f{t) = fc + fos{t)- For fos{t) 
we chose a square- wave signal of amplitude. A, and fre- 
quency, I'os = 1/Tosi where Tqs is the oscillation period 
(Fig. Id, upper). The four distinct levels of extension 
observed (Fig. Id, middle) correspond to the molecular 
extensions of the hairpin in the F and U states at the 
two force values, f = fc + A and f — fc — A. The power 
spectral density, S{i'), is defined as the Fourier transform 
of the stationary correlation function of the signal X{t): 



dt {X{t)X{0))c 



-i27ri>t 



(1) 



where (•) denotes a time- average over the signal. As 
shown in Fig. Id (lower), S{i') can be described as 
the superposition of a background power spectral den- 
sity, S'Ar(z^), and a structure of delta spikes centered at 
i^n ~ (2n-|- l)i/os ('^ = 0, 1, 2, • • ■ ). In order to extract the 
signal from the background noise, we define the output 
signal (OS), the background noise (BN) and the signal- 
to-noise ratio (SNR) as [Tl| . 



0S= lim 

Ai/->0 



BN = Sn{i^o 
OS 
BN 



SNR: 



S{v)dv. 



lim 



!>os + Aiy 



(2) 
(3) 

S{v)dv. (4) 



0.66 ± 0.04 s^i (Table SO in 



The SNR defined in Eq.Q is equal to the ratio of the 
spectral power of the signal at the frequency v^s (OS), to 
the noise-floor spectral density measured in the presence 
of the oscillation (BN) and has dimensions of Hz. Fig. 
Id (lower) illustrates how we measured the OS (red area) 
and the BN (blue vertical bar) from the s pec tral density. 
Other equivalent definitions of the SNR [17[ are the di- 
mcnsionless ratio between the power in the output signal 
(Eq.([2)) and the total input noise power delivered by the 
noise (proportional to the integral of background spec- 
tral density Sm{v) over all v). Because the total input 
noise power only depends weakly on v^s , we can take the 
OS, Eq.(I2]), as another indicator of the SR phenomenon. 
Indeed, both indicators OS and SNR are equally valid to 
identify resonant behavior, even thou gh t he peak is often 
more visible in the latter (see below) [T^ . 

For the hairpin HI at high trap power and trap stiff- 
ness Ktrap — 70 pN/^m, the resulting OS and BN as 
a function of i^os are depicted in Fig. 2a (lower), while 
Fig. 2c shows the SNR. In contrast to the OS, the pres- 
ence of a peak around Vcs ~ 0.4 ± 0.05 Hz is apparent 
for the SNR. This value is close to that predicted by the 
matching condition, z^sr = fcc/2, which states that the 
SNR is maximum when the average hopping time of the 
hairpin {1/kc = 1.56 s) is equal to half the period of 
the forcing oscillation (1/2j/os=1.25 s) [ll|, [ll-il ■ This 
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FIG. 2: SR experiments for hairpin HI at different 
trap stiffness. Pulling cycle (unfolding, blue; refolding, red), 
hopping trace, OS and BN for HI with amplitude A = 0.7 pN 
at high trap stiffness, Ktmp = 70 pN//im (a), and at low trap 
stiffness, Ktrap = 24 pN//im (b). Results averaged over 5- 
10 molecules. (Note: the force rips shown in force-distance 
curves should drop vertically without any finite stiffness cor- 
rection. The finite slope correction shown in Fig. 2b (top, 
left) is due to low-bandwith filtering of data), (c) SNR at high 
trap stiffness (low trap power) depicted in black (red). Units: 
OS (nm^), BN (nm^/Hz), SNR (Hz). Simulation results are 
shown as dashed lines (Figs. S4 and S7 in SI). The error bars 
represent the standard error over different molecules. 



shows that SR in singlc-moleculc hopping experiments 
approximately fulfills the matching condition as has been 
observed in other bistable systems. 

The OS and the SNR can be calculated theoretically 
as a function of the oscillation frequency for a Brownian 
particle in a continuous double- well potential [3 [3 113 ■ 
In this model, the OS and the SNR exhibit a soft and 
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FIG. 3: Molecular variability of the measured re- 
sponses (HI). Resuhs for the OS, BN and SNR for 10 

different molecules at high trap power Ktrap = 70 pN/^m. 
Units: OS (nm^), BN (nmVHz), SNR (Hz). The error bars 
represent the standard error over different molecules. 



a sharp peak, respectively, only when SR is induced at 
large enough forcing amplitudes [l^. These large forc- 
ing amplitudes correspond to a non-linear regime of the 
system, in which the shape of the double-well potential 
is so deformed that the barrier separating the wells van- 
ishes at the maximum elongation of the oscillation. In 
our experiments, we applied a large oscillation amplitude 
(^4=0. 7 pN). Note that the region of coexistence between 
the F and U states spans less than 3 pN in Fig. la (lower 
right). Thus an extra-force of 0.7 pN strongly alters the 
barrier and the relative free energy between states F and 
U. Our experimental results agree with the theoretical 
predictions by Stocks [l8| obtained in the non-linear re- 
sponse regime. We performed a numerical simulation of 
an overdamped particle moving in a double-well poten- 
tial with parameters that fit the experimentally measured 
molecular free energy landscape (Sec. IV in SI). Despite 
its simplicity, the model qualitatively reproduced the ex- 
perimental results for the OS, BN and SNR (dashed lines 
in Fig. 2c). 

In order to see what happens for lower oscillation am- 
plitudes, we explored the response of hairpin HI to an os- 
cillating force of lower amplitude, A = 0.2 pN. A very soft 
peak and a gentle maximum in the OS and the SNR can 
be seen around 0.4 Hz (Fig. SI in SI) in agreement with 
the results previously obtained for the higher amplitude, 
A = 0.7 pN (Fig. 2). However, the peak for A = 0.2 pN 
is much less clear than the peak for A = 0.7 pN, showing 
the importance of using oscillation amplitudes beyond 
the linear-response regime {AX^'^ /ksT ^ 1, where X''"'' 
is the characteristic distance separating the folded or un- 
folded states from the transition state. See also Sec. HI 
in SI for SR behavior in the linear response regime). 

A characteristic feature of SR experiments at the 
single-molecule level is the large variability observed in 
the measured response from different molecules. Fig. 3 
shows the OS, BN and SNR for 10 different molecules. 
Larger variability is observed for the OS as compared to 
the BN. This might be due to non-linear effects which are 
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sensitive to small differences in the experimental setup 
(e.g. tether misalignment with respect to the pulling di- 
rection, variations in the size of the bead and the trap 
stiffness, etc.). 



IV. INFLUENCE OF TRAP STIFFNESS AND 
LENGTH OF THE HANDLES 

An important issue in single-molecule experiments con- 
cerns the influence of transducing effects induced by the 
experimental setup (e.g. trap stiffness and length of 
the handles) on the measured kinetics. Recent studies 
[3 Ell, [13 show that the kinetic rates are only mod- 
erately affected (within one order of magnitude) when 
changing the length of the handles one thousand-fold or 
the trap stiffness ten-fold. Besides, numerical simulations 
carried out in Ref. show that kinetic rates for hair- 
pins measured with handles and trap always remain close 
and converge to the intrinsic rate (i.e. the rate measured 
without handles and trap) in the limit of very compli- 
ant linkers. To inquire the influence of the experimen- 
tal design on the kinetics of hairpin HI, SR was further 
investigated by varying conditions of the experimental 
setup such as 1) the stiffness of the optical trap and 2) 
the length of the handles. We observed how both effects 
changed the intrinsic noise of the system (Figs. 2b, 2c 
and 4). In the first case, when the trap stiffness, Ktrap, 
was decreased from 70 pN//im to 24 pN//im (Fig. 2b), 
the maximum peak in the SNR was shifted to higher fre- 
quencies (from 0.4 Hz to ~ 0.8 Hz) and became less clear 
(Fig. 2c, red curve). The effect of the trap stiffness on 
SR was evaluated by using the numerical simulation (Sec. 
IV in SI), finding good agreement between experiments 
and simulations (Figs. 2b and 2c). 

In the second case, if we increase by twenty-fold the 
length of the handles (528 bp and 874 bp at each flank- 
ing side) keeping the trap stiffness constant, Ktrap = 
70 pN/^m, we find that the resonance frequency shifts 
to a larger value for the long handles (Fig. 4). For the 
long handle construct, the matching condition is verified 
(i^SK = 2 Hz) and fcc ~ 4 s^^ as obtained from hopping 
experiments p^ . 

The dependence of the resonance frequency measured 
from SR, vsR, on the trap stiffness and the length of 
the handles was similar to that reported for the hop- 
ping rate measured in the hopping experiments at the 
co-existence force [isl . [isl ITg} . In both cases, the soft 
trap stiffness or the larger compliance of the long han- 
dles contributes to increase the hopping rate, supporting 
the conclusions of Rcf. [l^. Interestingly enough, the 
quality of the resonant peak worsens as the trap stiffness 
decreases but not as the linker becomes softer, show- 
ing that the quality of the SR peak is only dependent 
on the combined effective stiffness of bead and handles 
{k~^ = K^^l^ + K~l^^i^ ~ Kt~rip) which is approximately 
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FIG. 4: SR experiments for hairpin HI with long 
DNA handles, (a) Pulling cycle (unfolding, blue; refolding, 
red), hopping trace, OS and BN depicted in violet, (b) The 
resulting SNR in the case of high trap stiffness 70 
pN//im and the amplitude A = 0.7 pN. Results averaged over 
5 molecules. Units: OS (nm^), BN (nm^Hz), SNR (Hz). 
The error bars represent the standard errors over different 
molecules. 



V. OTHER SR QUANTIFIERS 

Next we investigated other representative SR quanti- 
fiers. These are: the fraction P\ of transitions that occur 
every half-period of the oscillation d, [H, [23| ; and the av- 
erage dissipated work, W [1, . To extract P\ , we mea- 
sured the residence time distributions, f (tf) and P(ru), 
of the F and U states in the presence of the oscillating 
force. The distributions are shown in Fig. 5a for hairpin 
HI in the cases i^os ~ 0.4 Hz (upper) and I'os = 5 Hz 
(lower) with yl=0.7 pN. Unlike the distributions shown 
in Fig. Ic, P(tf) {P(j\])) is not monotonically decreasing 
with rp (tu) and exhibits spikes corresponding to higher 
harmonics for rp = Ta^{\ -f 2n)/2 (ru = To^{\ 2n)/2) 
where n = 0,l,2,---. A few harmonic frequencies are 
shown as vertical arrows in Fig. 5a. In particular when 
i^os is close to the resonance frequency, the shape of the 
residence time distribution strongly deviates from an ex- 
ponential and a broad peak appears around the funda- 
mental mode, Tp = T'os/2 (ru = ros/2) (Fig. 5a, top). In 
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FIG. 5: Other SR quantifiers in hairpin HI. (a) Residence time distributions for the F (red) and U (blue) states at 
Vos = 0.4 Hz (upper) and Uos = 5 Hz (lower). The continuous black lines are the results of the simulations (Sec. IV in SI), (b) 
Pi in the F (red) and U (blue) and average work W (black) as a function of input frequency. Simulation results are shown by 
the dashed lines (Sec. IV in SI). Vertical dotted line shows the expected resonance frequency. The maximum in Pi is broad 
around 0.3-0.6 Hz, while that in W is also broad with a maximum found at a lower frequency in the range 0.1-0.4 Hz. (c) 
Average residence times (upper panel) and Pi (lower panel) in the F (red) and U (blue) states as a function of input frequency. 
Note that the frequency range is larger than that shown in b). Simulation results are shown as dashed continuous lines (Sec. IV 
in SI). The vertical lines show the two frequencies characteristic of stochastic resonance (vsr) and resonant activation {vra)- 
Statistics: 8 molecules, 3-5 minutes traces and 300-600 hopping transitions at each input frequency. The error bars represent 
the standard errors over different molecules. 



contrast, many peaks appear in the residence time dis- 
tribution when i/qs ^ i^SR (Fig. 5a, lower). 

Pi can be extracted from the area of the residence time 
distribution around the peak located at the fundamental 
mode, TF = ros/2 (ru = Tos/2). Let {n;i = 1, ... ,7V} 
be the series of N residence times measured in the pres- 
ence of the oscillating force. By counting the number, 
n, of Ti that satisfy the condition Tqs/S — Tos/4: < Ti < 
Tos/2 + ros/4, we define 



(5) 



Pi takes a large value if the residence time of the hairpin 
is equal to half the period of the oscillating force. This 
means that a large fraction of hopping transitions occur 
when the oscillating force changes sign. Therefore, the 
value of Pi has a maximum when SR is induced, because 
the transitions between the two states arc synchronized 
with the oscillating force (Pi is a bona fide SR quantifier 



[23| . See also Sec. Ill in SI). The results obtained for 
Pi in hairpin HI are shown in Fig. 5b. Pi exhibits 
a broad maximum around the resonance value z/gR. = 
kc/2 = 0.4 Hz. The broadness of the peak is in contrast 
to the narrower peak observed in the SNR (Fig. 2c). 
These results are consistent with analytical calculations 

For the average cyclic work done by an oscillating force, 
we define I26l 



w^-{jxdf) = {jfdX), 



(6) 



where the brackets stand for statistical averages over 
traces. Because W takes a large value when the fold- 
ing/unfolding of the hairpin is synchronized with the os- 
cillating force, it is a useful SR quantifier as well d, H^] . 
In fact, the larger the synchronisation between transi- 
tions of the hairpin and oscillations in the force, the larger 
the work done by the optical trap on the molecule. Re- 
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suits for W are shown in Fig. 5b. In contrast to SNR but 
similarly to Pi , the maximum in W is broad. Finally, we 
compared our experimental results with the predictions 
obtained from the numerical simulations in the continu- 
ous double-well potential whose parameters are the same 
as those used in Fig. 2 (Sec. IV in SI). Figs. 5a and 5b 
show a good agreement between experiments and simu- 
lations. Although both Pi and W show broad maxima as 
a function of Uos, they are not coincident, the maximum 
for the work is found at a lower frequency as compared 
to Pi. As pointed out in the introduction, the precise 
value of the resonance frequency depends on the quan- 
tifier specially when the quality of the resonant peak is 
low. 



VI. RESONANT ACTIVATION 

In stochastic systems driven by oscillating forces, it 
is customary to distinguish two effects: stochastic res- 
onance (SR) and resonant activation (RA). SR stands 
for the optimization of the response of the system (i.e. 
the output signal) whereas RA stands for the optimiza- 
tion of kinetics (i.e., maximization of the number of hop- 
ping transitions per second). SR and RA are different 
phenomena related to barrier crossing dynamics along 
temporally modulated energy landscapes RA is in- 
duced when the mean residence times of the states of the 
system are minimized with respect to the frequency of 
the oscillating force, i^ra- The values of vsR a-nd vha 
are often not the same, the latter being typically larger 
than the former. Fig. 5c (top) shows the mean residence 
times, (tp) and (tu), for hairpin HI measured in the 
range 0.1 Hz < t'os < 5 Hz. Only at higher frequencies 
(between 1 Hz and 2 Hz), the graph suggests a very shal- 
low minimum for the residence times. Therefore we are 
capable of observing both the SR and RA phenomena in 
the single-molecule experiments. The experimental re- 
sults also agree with the numerical simulations (Fig. 5c, 
dashed lines). Similar behavior has been reported in the 
experiments with a colloidal particle in a double- well po- 
tential generated by optical tweezers Q • 



VII. SR IN SHORTER HAIRPINS 

SR might be used to detect the transitions in cases 
where the hoppings of a hairpin is hard to be discrim- 
inated. These correspond to cases in which the hop- 
ping signal (extension jumps) are on the same order 
of the standard deviation of noise fluctuations. To in- 
vestigate this problem, we designed two short hairpins 
(SHIO and SH8) having only 10 and 8 base pairs along 
the stem, respectively (sequences shown in Figs. 6c and 
6d). The molecular free energy landscapes were calcu- 
lated for the two sequences at the theoretically predicted 
co-existence forces using the nearest-neighbour model for 
DNA (Fig. 6a, upper left) [H,!!^. As the length of the 



stem decreases, the landscapes show progressively lower 
co-existence force values, molecular extensions and ki- 
netic barriers. Measurements for SHIO and SH8 were 
taken at low trap stiffnesses to decrease the hopping sig- 
nal (/vtrap — 32 pN//xm and 17 pN///m, respectively). 
Pulling curves and hopping traces in the CFM are also 
shown in Fig. 6a (lower left). While the transitions are 
still visible for SHIO, these are hardly discriminated for 
SH8. This is also apparent from the dwell distributions 
on trap position, X, shown in Fig. 6a (right). Measured 
jumps in the molecular extension upon unfolding/folding 
are equal to 10.5 ±0.5 nm and 7.0 ±0.5 nm for SHIO and 
SH8, respectively. 

Fig. 6b shows the power spectra of X{t). Whereas 
SHIO can be fit reasonably well to a sum of two 
Lorentzians with two characteristic corner frequen- 
cies (0.64±0.02 Hz and 2.4±0.3 kHz), the quality of 
the fit considerably worsens for SH8 (~ 9.8 Hz and 
~ 15.6 kHz). The low frequency (in the range of Hz) 
in the power spectra corresponds to the hopping kinetics 
of the hairpin whereas the high frequency (in the range 
of kHz) corresponds to the random motion of the opti- 
cal trap caused by the force feedback. Because the noise 
in the trap position, X, introduced by the force feed- 
back protocol is not of thermal origin, the power spectra 
measured in the CFM should not necessarily be fit to a 
sum of two Lorentzians. This is specially acute for SH8 
where the feedback loop cannot follow the fast hopping 
transitions. 

Once the hopping properties of the hairpins were char- 
acterized, we then carried out the oscillating experiments 
for hairpins SHIO and SH8 around the co-existence force. 
The results we obtained for SHIO are similar to those re- 
ported for hairpin HI at low trap power shown in Fig. 
2c. For SHIO the peak in the SNR around i'sr=0.5 Hz is 
close to fcc/2 where kc was measured to be 0.43±0.07s 
from the hopping traces for X{t). More interesting is 
the case of hairpin SH8 where the co-cxistence force can 
still be located, but the hopping signal is blurred by the 
fluctuations. In Fig. 6d, we can see that the OS and the 
SNR exhibit a maximum around vsR = 5 ± 1 Hz for SH8 
which gives fee ~ 10 ± 2 s"""^ according to the matching 
condition. This value agrees with the value of ~ 9.8 Hz 
obtained from the Lorentzian fit to the power spectrum. 
As an additional test, we have implemented a Hidden 
Markov Model (HMM) with the forward-backward feed- 
back algorithm as described in Ref. [lOl to extract the 
kinetic rates of SH8 from the hopping trace, X{t). By 
applying the HMM to the hopping traces of SH8, we ob- 
tained a value of kc = 9.4 ± 0.5 s~^ (7 molecules), which 
confirms the results obtained with SR and Lorentzian fit 
to the spectral density. 

Thus, SR offers an alternative method to estimate the 
hopping rate of SH8. Indeed, the two states (F and U) 
cannot be easily detected from the hopping trace and the 
residence time analysis done for hairpin HI (Fig. Ic) is 
difficult to implement. In this case SR confirms the value 
of the hopping frequency initially obtained from a poor 
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FIG. 6: SR experiments in shorter hairpins, (a) Free energy landscapes (upper left), force-distance curves (lower left) 
and hopping traces in the CFM for SHIO (red), SH8 (blue) and HI (black). Measurements were carried out with a low trap 
power Ktrap — 32 pN/fim for SHIO, 17 pN//im for SH8, respectively. For sake of simplicity, all pulling curves in the lower 
left panel are shown parallel with equal average slopes. The dashed-squared region for SH8 curves indicates the region where 
unfolding/folding transitions occur. Distributions of trap position, X, show clear transitions for SHIO but not for SH8. (b) 
Power spectrum for SHIO and SH8. Cyan curves are fits to a sum of two Lorentzians (see text for details). Colors as in Fig. 
6a. (c, d) OS, BN and SNR for SHIO and SH8. The amplitudes of oscillation force are A = 0.5 pN for SHIO and A = 0.15 pN 
for SH8, respectively. Colors as in Fig. 6a. Statistics (SHIO, SH8): Molecules (5, 7); Duration of hopping traces (4, 2) minutes; 
Average number of hopping transitions (250, 1200) at each input frequency. Units: OS (nm'^), BN (nm^/Hz), SNR (Hz). The 
error bars represent the standard errors over different molecules. 



Lorcntzian fit of the power spectrum. 



VIII. CONCLUSION 

We carried out SR experiments in single DNA hair- 
pins subject to an oscillatory mechanical force of varying 
frequency. Our aim was to investigate how a molecule 
exhibiting bistability (i.e. hopping between the folded 
and unfolded conformations) responds to an applied os- 
cillating force. In SR the response gets amplified at fre- 
quencies close to the characteristic hopping frequency of 
the hairpin. By measuring the power spectral density 
of the molecular extension, we carried out a detailed in- 



vestigation of the frequency dependence of the output 
signal (OS, Eq. (2)), the background noise (BN, Eq. 
(3)) and the signal-to-noise ratio (SNR, Eq. (4)) in the 
20bp hairpin HI which exhibits dichotomous hopping 
behavior. We then extended our research by exploring 
how several parameters of the experimental setup such 
as trap stiffness, length of the handles, oscillating am- 
plitude and size of the hairpin influence the resonance 
behavior. From the measured traces, we also analyzed a 
few other SR quantifiers such as the number of folding 
and unfolding transitions occurring every half-period of 
the oscillation (Pi, Eq. (5)), the average mechanical work 
per period of the oscillation (W, Eq. (6)) and the mean 
residence times in the unfolded and folded states ((tu) 
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Comparison between usr (Hz) and fee (s ^) 





t-'SR iiOlli oiNlX 


i-'SR from OS 


i^SR from _Pi 


t^SR from 


i- /9 


HI"'" 


0.40±0.02 (n = 10) 


0.45±0.03 (n = 10) 


0.45±0.04^ 0.48±0.08' (n = 8) 


0.33±0.07 (n = 8) 


0.33±0.02 (n = 12) 




0.72±0.08 (n = 5) 


0.50±0.06 (n = 5) 






0.53±0.07 (n = 8) 




2.0±0.2 (n = 5) 


1.9±0.1 {n = 5) 






2.2±0.3 (n = 5) 


SHIO 


0.54±0.02 (n = 5) 


0.58±0.04 (n = 5) 






0.43±Q.07 (ji = 4) 


SH8 


5.7±0.4 (n = 7) 


5.3±0.2 (n = 7) 






4.7±0.3'=^ (n = 7) 



TABLE I: Resonance frequency, vsn, obtained from SNR, OS, Pi and W vs. flopping rate, kc, at the co-existence force (see 
Sec. I in SI for kc). i^sk was ciiosen as tiie peak value of eacii SR quantifier for eacli molecule, (a. Short handles, b. Long 
handles, c. High power trap (k — 70 pN//im). d. Low power trap (k — 24 pN//im). e. Folded state, f. Unfolded state, g. 
Rate determined using a Hidden Markov Model (HMM). Note that n is the number of molecules analyzed.) 



and (tf)). The mean residence times describe a mech- 
anism slightly different from SR that has been termed 
resonant activation (RA). Overall, we find that the SNR 
and the other SR quantifiers (such as OS, Pi, W) ex- 
hibit a peak at a frequency close to that determined by 
the resonance matching condition. Among all quantifiers 
only the SNR and the OS tend to show a modest amplifi- 
cation of the response, the SNR showing a higher quality 
peak. Our results are summarized in Table L Moreover, 
our experimental results are well predicted by numerical 
simulations of an overdamped particle in a double-well 
potential reproducing the measured molecular free en- 
ergy landscape of the hairpin (Sec. IV in SI). Finally, 
our experimental findings also agree with theoretical re- 
sults [l8| that show a modest gain in the response of noisy 
systems driven by oscillating forces. 

A unique aspect of our work is the investigation of SR 
in small systems in conditions of weak thermodynamic 
stability (folding free energies of a few k^T units) not far 
from noise level (fceT). This has a primary consequence: 
the proper control parameter in our experiments does 
not appear to be the noise intensity. In fact, by chang- 
ing noise intensity (e.g. by tuning temperature or dc- 
naturant concentration), we also modify the structural 
properties of the molecule in a non-controlled way (i.e. 
by changing its thermodynamic stability or free energy 
of formation). Our work circumvents this problem by 
using the frequency of the external driving force as con- 
trol parameter. Simple as this choice may seem only a 
few theoretical and experimental works have addressed 
it in the past. From this perspective, our study should 
stimulate further theoretical work in SR of small sys- 
tems where noise intensity and thermodynamic stability 
are tightly coupled. Another consequence of the noise in- 
tensity vs thermodynamic stability coupling is the strong 
variability exhibited by single-molecule SR experiments: 
the measured signal-to-noise ratio versus any control pa- 
rameter (in our case, oscillation frequency) will tend to 
show large variations from molecule to molecule. This 
was apparent in the results for hairpin HI shown in Fig. 
3 and has been observed in the rest of molecules (see, for 
instance the results shown in Fig. 7 for SH8). Such vari- 
ability is consequence of the aforementioned weak stabil- 
ity of biomolecular bonds, and various sources of exper- 
imental errors (e.g. instrumental drift, misalignment at- 
tachment, inaccurate discrimination of the co-existence 



force, etc.). It has no counterpart in other SR stud- 
ies of non-linear macroscopic devices or single degree-of- 
freedom systems (such as single colloidal particle in op- 
tical traps or macroscopic systems in solid state physics 
or electronic devices). 

IX. FUTURE PERSPECTIVES 

The results of our work suggest that we could extract 
the kinetic rates of molecular hoppers by measuring the 
resonance frequency in oscillating experiments. Is this 
approach useful? There are several widely accepted and 
commonly used single-molecule methods that can ex- 
tract the kinetic parameters of molecular hoppers just 
by analyzing the hopping traces without bothering about 
carrying out oscillating measurements. It is then clear 
that single-molecule SR is not worth pursuing if other 
simpler methods are available. Yet SR might be of in- 
terest for investigating fast molecular transitions where 
current methods might fail. In Section VII, we investi- 
gated SR in an 8bp short DNA hairpin (SII8) at condi- 
tions (low trap stiffness) where hopping rates are hard 
to be measured from standard methods (e.g. the Bell- 
Evans model). The faster hopping rate and the smaller 
jumps in extension (due to both the shorter length of 
SH8 and the decreased trap stiffness) contribute to make 
the hopping rate measurements difficult. Note that we 
have been able to extract the value of the hopping rate 
either by measuring the power spectrum (Fig. 6b) or 
by implementing a hidden-Markov model. Interestingly, 
whereas applying standard methods to extract kinetic 
rates become steadily difficult as the hopping signal be- 
comes more noisy, the quality of the resonant peak in the 
SNR remains acceptable (Fig. 6d). This suggests that in 
experimental conditions where hopping signals become 
nearly undetectable, SR may find a fertile ground for 
useful applications. 

Measuring the kinetics of single bonds might be cru- 
cial to dissect the kinetic pathways of many reactions, 
from nucleic acid translocases indispensable in virtually 
all tasks of nucleic acids metabolism, to molecular folding 
of proteins and ligand-receptor binding. Moreover, the 
detection of single bond kinetics also provides a direct 
measurement of the affinity (or free energy of formation) 
of weak single bonds (e.g. important for an accurate 
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FIG. 7: Molecular variability of the measured re- 
sponses (SH8). Results for the OS, BN and SNR for 7 
different SH8 molecules. Units: OS (nm^), BN (nm^/Hz), 
SNR (Hz). The error bars represent the standard error over 
different molecules. 



determination of the parameters characterizing the ther- 
modynamics of secondary structure formation in nucleic 
acids [3^). It is therefore important to explore new ap- 
proaches capable of illuminating into such questions. The 
experimental resolution of formation/dissociation kinet- 
ics is currently limited to 5 base pairs [2^ [sij. Over- 
coming this limit strongly relies not only on increasing 
the hopping signal relative to the noise but also in slow- 
ing down the (expected too fast) formation/dissociation 
kinetics of single bonds. A direct measurement of the 
formation/dissociation fast kinetics of single molecular 
bonds stretchable along sub-nanometer scales and resis- 
tant to low (a few piconcwton) forces remains an exper- 
imental challenge. 

In fact, the route to discriminate hopping kinetics in a 
small number of base pairs may be plagued of difficulties. 
The situation might be even worse if the aim is to detect 
the unraveling kinetics of a single nearest-neighbor base 
pair (NNBP), which is the minimal unit of DNA bonds 
(double stranded helices are stabilized by both hydro- 
gen bonds between complementary bases and stacking 
between NNBP). Currently most kinetics measurements 
are carried out in hopping experiments. However there 
is a complication present in hopping experiments due to 
the low signal-to-noise ratio inherent to unraveling a sin- 
gle NNBP together with the disturbances caused by the 
multifrequential noise present in the high frequency range 
where the kinetic rate of formation/dissociation of a sin- 
gle NNBP is expected to fall. The low signal-to-noise 
ratio problem can be partially resolved using advanced 
data analysis tools such as Bayesian methods and HMM 
to unravel hopping traces for SII8. However such meth- 
ods assume a specific form of the noise (i.e. decorre- 
lated force fiuctuations and Gaussian emission signal) 
and do not account for multifrequential sources of noise 
(due to the aforementioned sources). In this regard, SR 
might be extremely useful to separate the true forma- 
tion/dissociation kinetics of a single NNBP from these 
other artifacts. 

Finally our work focused in the SR phenomenon 



in DNA hairpins wehereas other interesting molecular 
structures are now available for single molecule pulling. 
From this point of view, it would be very interesting to 
carry out SR measurements in more complex molecular 
folders (e.g. exhibiting multiple folding pathways, 
intermediate states or non-cooperative transitions) such 
as RNAs and proteins. 

Methods. 

Synthesis of DNA hairpins. The DNA hairpins with 
handles are synthesized using the hybridization of three 
different oligonucleotides (Fig. la). One oligonucleotide 
contains the sequence of the ssDNA left handle plus 
a part of the sequence of the desired DNA hairpin; 
the second has the rest of the sequence of the DNA 
hairpin and the ssDNA right handle. The right and 
the left handles have the same sequence in order to 
hybridize them with the third oligonucleotide. The first 
oligonucleotide has a biotin at its 5' end and the second 
oligonucleotide has been modified at its 3' end with 
a digoxigenin tail (DIG Oligonucleotide Tailing Kit, 
2nd generation, Roche Applied Science). Once the first 
and the second oligonucleotides are hybridized to form 
the hairpin, the third oligonucleotide is hybridized to 
the handles to form the dsDNA handles. Streptavidin- 
coated polystyrene microspheres (1.87 ^m; Spherotech, 
Livertyville, IL) and protein G microspheres (3.0-3.4 /im; 
G. Kisker Gbr, Products for Biotechnology) coated with 
anti-digoxigcnin polyclonal antibodies (Roche Applied 
Science) were used for specific attachments to the DNA 
molecular constructions described above. Attachment 
to the anti-digoxigenin microspheres was achieved first 
by incubating the beads with the tether DNA. The 
second attachment was achieved in the fluidics chamber 
and was accomplished by bringing a trapped anti- 
digoxigenin and streptavidin microspheres close to each 
other. The sequences of the short hairpins are: SHIO 
(5'-GCGGCGCCAGTTTTTTTTCTGGCGCCGC-3'), 
SH8 (5'-GGCGCCAGTTTTTTTTCTGGCGCC-3'). 
Experimental setup. The experiments have been 
carried out using a high stability newly designed minia- 
turized dual-beam optical tweezers apparatus j3^ . It 
consists of two counter-propagating laser beams of 
845 nm wavelength that form a single optical trap where 
particles can be trapped by gradient forces. The DNA 
hairpin is tethered between two beads (Fig. la). One 
bead is immobilized at the tip of a micropipettc that is 
glued to the fluidics chamber; the optical trap captures 
the other bead. The light deflected by the bead is 
collected by two photodetectors located at opposite sides 
of the chamber. They directly measure the total change 
in light momentum which is equal to the net force acting 
on the bead. Piezo actuators bend the optical fibers and 
allow the user to move the optical trap. The force is 
made to oscillate using a force feedback system that op- 
erates at 4 kHz minimizing instrumental drift effects as 
compared to protocols without feedback. Force feedback 
does not introduce artifacts in our measurements unless 
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i^os is too high (typicahy larger than 50 Hz) or A is too 
small (less than 0.1 pN). 

The folding-unfolding experiments described in this 
report were performed at room temperature (24°C) in a 
buffer containing lOmM Tris-HCl pH7.5, ImM EDTA, 
IM NaCl, 0.01% Sodium Azide. 
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